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1 Introduction 


Let X = (X(t): t £ [0, T]} be a second order stochastic process with the incremental variance 
function o\ defined on [0,T ] 2 := [0, T\ x [0, T] with values 

a: x (s,t) :=E[X(t)-X( S )] 2 , (s,t) € [0,T] 2 . 


Denote by $ a class of continuous functions ip: (0,T] —» [0, oo) such that lim/qo+(h) = 0 
and L(h) = ip(h)/h —¥ oo, h J. 0. For example, we can take <p{h) = h\ lnh| or p(h) = h 1_e for 
0 < s < 1. Set 


7* : 


7* = 


f h' 1 

= inf < 7 > 0 : lim sup --— 

l ip(h)^s^T-h CX[S, S + h) 

( h 7 1 

inf \ 7 > 0 : lim-—— = 0 1 

l hro (j y (0, h) J 


= 0 


and 


:= sup < 7 > 0 : lim 
[ h.10 

:= sup < 7 > 0 : lim 

l hlO 


inf 


h? 


hiO ip(h)^s^T-h ax(s, s + h) 

■}. 


hio cr x (0, h) 


= +00 



where tp G $. Note that 0 ^ 7 * ^ 7 * ^ +00 and 0 ^ 7 * ^ 7 * ^ + 00 . In paper [13] we 
used a narrower class of functions *&, i.e. such functions which additionally satisfy condition 
lim/^op • L 3 (h)] = 0. This condition is not necessary for the existence of the Orey index. 
It is required only for consideration of almost sure asymptotic behavior of the second-order 
quadratic variations of Gaussian processes. 

We give the following extension of the Orey index. 
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Definition 1 f [13| 1 Let X = (A'(t): t G [0,T]} be a second order stochastic process with 
the incremental variance function a\ such that sup 0<s<T _ h ax{s,s + h) —> 0 as h —> 0. If 
7 * = 7 * = 7 * = 7 * for any function ip 6 f, then we say that the process X has the Orey index 
Jx = 7* = 7* = 7* = 7* • 

Assume that for some 7 G (0,1) the second order stochastic process X satisfies conditions: 
(Cl) < 7 .y(0, 5) x 5 1 , i.e., ax{ 0, 5) and S 1 are of the same order as <5 | 0; 

(C2) there exist a constant x > 0 such that 


m 


sup sup 

ip(S)^t^T-5 0 <h^S 


ax(t,t + h) _ 
xh~< 


as S 4 . 0 


for every function <p G 'h. 

If for some constant 7 G (0,1) the second order stochastic process X satisfies conditions 
(Cl) and (C2), then the Orey index is equal to 7 (see |13|). 

Recently much attention has been given to studies and applications of Gaussian processes 
such as fractional Brownian motion (fBm), sub-fractional Brownian motion (subfBm), bifrac- 
tional Brownian motion (bifBm), fractional Ornstein-Uhlenbeck process. All of them are 
Gaussian processes and they have the Orey indexes. Consequently, examining Gaussian pro¬ 
cesses with the Orey index we thus examine the processes listed above. 

Many authors (see m, hd], m, s, el m , ra, hzi, usd considered an almost sure 
convergence and asymptotic normality of the generalized quadratic variations associated to 
the filter a (see m of a wide class of processes with Gaussian increments. The strong 
consistency of the Orey index estimator was proven in ns. 

In the papers of Breton et al. [3] and Breton and Coeurjolly [5], an exact (non-asymptotic) 
confidence interval for the Hurst index of fBm was derived with the aid of concentration 
inequalities for quadratic forms of Gaussian process. The obtained confidence intervals for the 
Hurst parameter were based on a single observation of a discretized sample path of the interval 
[0,1] of fBm. Exact confidence intervals for sub-fractional Brownian motion were considered 
in m but are not sufficiently precise. 

The purpose of article is to extend the results of Breton et al. [3] and Breton and Coeurjolly 
[5] as well as to apply them to Gaussian processes with the Orey index which may not have 
stationary increments. 

The paper is organized in the following way. In Section 2 we give exact confidence intervals 
for the Orey index of Gaussian process. Section 3 contains some application results for known 
Gaussian processes which may not have stationary increments. Finally, in Section 4 some 
simulations are given in order to illustrate the obtained results. In addition, Appendix includes 
the R code listings of simulations. 


2 Confidence intervals 

First, we formulate a concentration inequality for a family of Gaussian r.v.’s. Consider a finite 
centered Gaussian family X = {X *. : k = 1,...,A/}, and write d^j = E (X^Xj). Define two 
quadratic forms associated with X and with some real coefficient c: 

M M 

Qi{c,X)=c^2(Xl-d k k), Q 2 (c,X) = 2 c 2 XkXjdkj- (1) 

k =1 fc,i=1 

The following statement characterizes the tail behavior of Qi(c, A'). 

Theorem 2 ( [2], S3 ) Suppose that Q\{c, A) is not a.s. zero and fix a > 0 and fi > 0. Assume 
that Q 2 (c, A) < qQi(c, A) + (3, a.s.-P. Then, for all z > 0, we have 

P(Qi(c, A) > z) ^p r (z;a,/3) := e _z/ “(l + ’ 

P(Qi(c, A) < -z) ^ipfiz; a,/3) := e z/a (l - —J 1 [a,p/ a ]{z). 

Remark 3 J^jj Note that ip r (-;a,/3) (resp. </?;(■; a, /3)) is a bijective function from (0,+ 00 ) 
(resp. (0 ,P/a)) to ( 0 , 1 ). 
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Next, we apply the obtained concentration inequality to second order quadratic variations. 
Let X be a a centered Gaussian process satisfying conditions (Cl) and (C2) with the Orey 
index 7 £ (0,1). Denote Y = ■ ■ ■, Y n ~i,„), where 


1 k.n 


T.xV = X{tl+l) ” 2X{tl) + 


t% = — and x is a constant defined in condition (C 2 ). Set d-jf 1 = E Yj, n Yk,n- 

Proposition 4 Assume that there exists a sequence of real numbers (e„) not depending on 7 
and such that 




< S n 


£n -i- 0 . 


( 2 ) 


Then for all z > 0 we have 


n — 1 \ 

E i( Y ^) 2 - 4r] ^ zj ^<fir,n(z;l'n,En), 

L-— 1 / 


( ~ d k ’k ] ^ — ^ ^<Pl,n(z',V n ,£ n ), 


where 

Tr,n(z, l^ni En) - = 6 Xp^ 
<Pl,n(z\ v n ,E n ) := exp 


— 1 


2u n 
— 1 

2 

n — 1 


1 + 


(^n + 1 )\/ 71—1 


( eTt + l)(n-l) 


1 - 


(£n + 1 ) n — 1 


(e„ + l)(n-l) 


L [0,(e Tl + l)\/n 


=!](*)> 


""=EKfi- 


t=i 


Proof. Denote 


and 


Qi((n -1 r 1/2 , y) = 2 [(n ,„) 2 - <&"] 

vn 1 fc=i 


Q 2 ((n - 1)- 1/2 ,P) = E ^ <?■ 


fc.3=l 


Then 


Q 2 ((n-1)- 1/2 ,P) 

n-t n-1 

<ittE i^-i • 1^1 ■ Krl < —r E K 2 + *?] • l d £ 


k,j = 1 
n — 1 


n — 1 

fc,j=i 

n—l n—1 


= E n 2 ME"i = E y * E \ d kf 


k,j = l 

n — 1 


k=1 j=1 


n —1 


2 z/ n 


n — 1 

E^U ^EKf 

i=i 


= ^((n - 1)- 1/2 , F) + ^=L= ]T [dlj? - 1] + 


y/n — 

^a„Qi((n-l)“ 1/2 ,y)+A 


where 


a„ = 2Vn , /3„ = 2v n {e n + 1 ). 

Vn — 1 
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(3) 

(4) 


Thus 

P(Ql((n- 1 )~ 1/2 ,Y) ^ z) ^ <Pr,n{z;Vn,£n) 

P(Qi((n- 1 )" 1/2 ,F) ^ -2) ^ 1 pi,n(z;u„,£ n ) 
and the proof is completed. ■ 

Remark 5 Since ipi, n (z; £ n ) and <p r ,n(z; -,e n ) are non-decreasing, inequalities 0 and 0 
remain true with v n replacing v n , where Vn^v n . 

For any a G (0,1) and v n (y n ^ v n ), denote by qi, n {a) := (y3i,„) _1 (a; v n , e n ) and q T ,„(a ) := 
(ip r ,n)~ 1 {or,Vn,£n)- For convenience we define 


Xl 


/ \ 1 QZ,7i —l(a) , / \ i Qr,n— l(a) 

i-i(a) := 1 - and *r,n-i(a) := H- , , ■ 

— 1 v n — 1 


Note that Remark[3]above ensures that for any a G (0,1) and for all n > 1, we have xi, n -i(a) > 
0. Set 

n — 1 

Sn := J2 ( A( Z X ) 2 ’ 9 u,t(x) := 2*In (n/T) - ln(4 - 2 2 *), * G (0,1), 


and 


In* x = 


—00 if x = 0, 

In a; if x > 0. 

The function gn,r(x) is a strictly increasing bijection from (0,1) to (— In 3, + 00 ) if n > T. 
Theorem 6 Let a G (0,1). Assume that conditions of Proposition^ are satisfied and there 
are constants v„ such that v n ^ v n . Then 


P(7 G [7n (a)>7n“ P («)]) > 1 - 01 , 


where 


7n n/ (a) := max ( 0 , g\ n ' T max ( In* ( ( — ^ V 0^ , - In 


7 u up (a) := min 1 ,g, 


9n,T 


In 


n-l («/2) + E n 2 

s * 

On 


Proof. Denote 

Z n = (n- l)- 1 / 2 R ?l ( 2 ) (y,2) - V^l = v^l[(n- l)- 1 V^ 2 ) (y,2) - 1], 

where 

n — 1 

K C2) (P,2) = ^(A^y) 2 . 

4=1 

Then 

| - q f Z,n-l(o;/2) < (n - 1) _1/2 ^ [(Yfc,n) 2 - d^. n ] ^ 5r,n-l(a/2)| 

^ fc=l ' 

{ n—1 n —1 n 

- n-i(a/2) + (n - 1)^ 1/2 ^ [dfei" - l] ^ < ?r,n-i(a/2) + (n - 1) _1/2 ^ [dfcfc" - l] !■ 

fc=l fc=l ' 

{ n-1 27 

xi,„~i{a/2) + (n - I)" 1 ^ [dj£* - 1] < y27 _ 47) S n 

n — 1 \ 

< Xr,n -1 (a/2) + (n - 1) _1 ^ [dl’ k n - l] l 

fe=i ' 


,27 


G s a;i,n-i(a/ 2 ) s n ^ 


T 2 7 >f 2 ( 4 _ 227) 
27 


S/i ^ X r ,n — l(a/2) -(- £7 


f n 2 7 'I 

= | (xi, n-1 (a/2) -£„)v0^ y 27><2 ( 4 _ 2 2 7 ) -S'n 71-1 (a/2) + £n j 

= { In* ((®(,„_i(a/2) - £„) V0) -lnS„ + lnxr 2 ^ g n , T ( 7 ) < ln(a; r ,„_i(a/2) + £ n ) - lnS 1 ™ + lnx: 2 


In* 


(xi,n-l(a/‘2) — £ n ) V 0 2 

-rr- yt, 


V (-In3) < 9ti,t( 7 ) < In 


%r,n -1 (^/2) + £n 2 

-7“- ^ 
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Note that 


p( - ®,„_i(a/2) ^ ~~~~ ^ [(n,n) 2 - dfci"] ^ ? r ,„-i(a/2) > ) > 1 - a. 

' V n 1 k=1 J 

Thus 

P(t£ [7r / («),7r p (a)0 ^ 1 - a. 

The proof is completed. 


3 Applications 

In this section we obtain the confidence intervals for subfBm, bifBm and the fractional 
Ornstein-Uhlenbeck process. For this purpose we apply the Theorem [ 6 ] In order to apply 
the Theorem [ 6 ] it suffices to find the sequence of real numbers (e„) in the estimation © and 
estimate v n . In the considered cases as the special case appears the Brownian motion. We ex¬ 
clude it from consideration in view of its properties (in particular, independent increaments). 
It is easy to see that 


^ 7 B ,71 I 

a kk ~ T 


n 

— max 
2 T 


l,£i4" 


< 2 , 


3 = 1 


where d 


' B , n — tpaI 2 ) ra< 2 ) 


'jk 


EAi /BA/ n B and B is Brownian motion. 


3.1 Sub-fractional Brownian motion 

Definition 7 ([3]) A sub-fractional Brownian motion (subfBm) with the index H , H £ 
(0,1), is a mean zero Gaussian stochastic process S H = (St 1 ,t ^ 0) with the covariance 
function 

G h (s, t) := s 2H + t 2H - i [(a + t) 2H + \s - t \ 2H ]. 

The case H = 1/2 corresponds to the Brownian motion. For H 1/2 this process has some 
of the main properties of fBm, but its increments are not stationary. 

The incremental variance function of subfBm is of the following form 

ct s h(s, t) = E|S t - S s \ =|t-s| +(s + t) -2 (t +s ). (5) 

For any 0 ^ s ^ t ^ T the inequalities (see [3]) 

(t-s) 2H sj a 2 sH (s,t) sj (2-2 2H ~ 1 )(t-s) 2H , if 0 < H < 1/2, ( 6 ) 

(2-2 2H ~ 1 )(t-s) 2H ^o 2 sH (s,t) ^ (t-s) 2H , if 1/2 < H < 1 (7) 


hold. 

It is known that for subfBm the Orey index is equal to H (see 031)- Now we prove the 
following lemma. 

Lemma 8 Assume that S H = {S H (t) : t £ [0,T]} is a subfBm. If H 7 ^ 1/2 then 


n — 1 


E d : 



In 33 \ 
6 T n- 1 + 9In4 ) ’ 


where 

„ J2H 

jS H ,n _ n _p,/ A (2) q H \ 2 

a kk — J~2H^ _ 2^ h ) V k,n J ) ■ 

Proof. Observe that the following equality 

E (S? +h - 2 S t H + S?_ h ) 2 =(4 - 2 2H )h 2H - 2 2H ~\t + h) 2H - 3 • 2 2H t 2H 

- 2 2H ~ 1 (t - h) 2H + 2(2t + h) 2H + 2(2 1 - h) 2H 
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holds. Thus 


,S H ,n _ 


1 - 


(k + l) 2 " + 3 • 2 2H k 2H + 2 2H ~ 1 (k - lf H - 2(2 k + 1) 2H - 2(2 k - l) 2 


4-2 2H 


_ b(k,H) 

4 - 2 2H ' 


oH 

For simplicity we shall omit the index n for ’ . Using computer modeling we obtain the 
inequalities 


7 \b(k,H)\ ^ \b(l,H)\ ^ 1 

max d kk ^ max dn = - and max -- — w ^ max -- —prr ^ x 

hs(o,i) hs( o,i) 6 ffe(o,i) 4 — 2 2iJ He(o,i) 4 — 2 2H 6 

Let ip £ 4' and denote r„ := [p(Tn~ 1 )nT~ 1 ], where [a] is an integer part of a. Then 


(8) 


n— l 

k- 

=h + h. 


\b(k,H)\ 


n — 1 


+ max 




n-1^4-2 2 ^ T„+i^fc<n-i 4 - 2 2ff 

k=1 


It is clear that 


r ^ r„ \b(k, H)\ ^ t„ 1 6(1, ff)| ^ 1 frr _ lw _! n 

A ^ -7 . max '/ „ 0 (/ -x x —^ x v{Tn )T 


n—% i 4-2 2H n- 1 4 - 2 2W 6 

Now we estimate 6(/c, FT) for k ^ 2. Using the formula 


n—l 


(i+*r = i+E 


a(a — 1) • • • (a — k + 1) k 
k\ * 


for — 1 < x < 1, 


we obtain 


b(k,H) =2 2 " _i (fc+ l) 2 " + 3 • 2 fc + 2 2M ~ 1 (k — l) 2 " - 2(2fc + l) 2 " - 2(2 k- l) 2 


=2 2H ~ 1 k 2H 


=2 2H ~ 1 k 2H 


1 


1 + t +6+11—— —4 1+— — 41 1 — — 


1 


E 


k 


2H(2H - 1) • • • (2ff - m + 1) / 1 


2fc 


2fc 


+ (-l) m - 2 2_m - (—l) m 2 2-m ] 


=2 2H k 2H 


^ 2H(2H - 

- 1) • • • (2If - 2 m + 1) / 

1 V m (i - 2 2 ~ 2m ) 

2 ^ 

- m =2 

(2m)! 

k) ^ ’ 


Note that the sign of 2H(2H — 1)(2 H — 2) • • • (2 H — 2m + 1) is the same as that of 2 H — 1. 
Thus 


\b(k, H)\ ^2 2H k 2H 
4-2 2H E- 2 2H 


E 


2H\2H — 11(2 — 2H) • • • (2m — 1 — 2H) (1\ 2 ™ U n2 _ 2ml 
(2m)! \k) [ [ 


3-2 2H k 2H ^ 2(2m — 1)! /1\ 2 ™ 


E 


4 In 4 ' (2m)! 

m=2 v 1 


©■ 


q c\2H 


4 In 4 
3- 2 2 


4 In 4 


v - 1/1 \ m 

•E m \k 2 ) ^ 

m =2 

(>-*) 


Q r\2H 

*5 ■ z 7 2H — 2 

~4hU k 


E l /I \™-i 

m — 1 V fc 2 / 


k - In (1 — +r I < 


33 • 2 2iJ 
36 In 4 




33 

9 In 4 


2ff|l - 2ff|(2 - 2H) ■ ■ ■ (2m - 1 - 2ff) s; 2(1 - if)(2m - 1)!, 


log(l — u) = — V' if 0 < w < 1 and — log(l — u) ^ if 0 < u ^ 1/4 

' At 9 


(1 — if) In 4 < 4 — 2 2H . 
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So 


- 


--lvn-1 n 


33 


^ 6 <p{Tn )T n _ 1 + gin 4 ( Tn + 1)2 

1 ! ! _33_ T 2 

^6 n—1 9In4 (^ 2 (Tn _1 )n 2 


Let ip(Tn 1 ) = {Tn 1 ) 2/ ' 3 . It belongs to the class of functions 4b After putting into the 
obtained inequality we get 


1 


n — 1 


1 




1 


6 T 1 / 3 n 2 / 3 n — 1 9 In 4 


33 /T\ 2/3 _ /T\ 2/ V 2 


2/3 


n 33 
+ 


n J \15T n- 1 91n4/ 


Lemma 9 Assume that S H = {Sp ■ t £ [0,T]} is a subfBm. If H Y 1/2 f/ien 

1^2’ 

3 = 1 


(9) 


where 


,s H ,n = n 2H jsH n 
a jk rp2H ^4 _ 2 2 ^) ’ 


= EA^A^S". 




Proof. The fourth order mixed partial derivative of the covariance function Gh{s, t) is of the 
following form 


8TG h 

ds 2 dt 2 


(s, t) = — Ch 


+ 


|s-f| 2 (2 — H ) (g + t^-H) 

for each s,t > 0 such that s ^ t, where Ch = H(2H — 1)(2 H — 2)(2 H — 3). Since the 
covariance function Grr(s,t) is continuous in [0, T] 2 and the derivative g %g t 2 is continuous in 
(0,T] 2 /{s = t} then for H Y 1/2 and j ^ 1 or k ^ 1 

~h r*?+i l~ u /•*£+1 r x 8 4 Gh 

d jk ’ = / du dv dx g^{v,y)dy. 

Jt 7 Ju-T/n Jt " Jx-T/n us ul 


Assume that Rh is the covariance function of the fBm B H . Then the derivative 

8 4 R H , .. C H 

r(M) = -7 


ds 2 dt 2K,J \s-t\ 4 ~ 2H 

of the covariance function Rh is continuous in (0,T] 2 /{s = t } and 
df k H ’ n :=E(A ™B h ^ 2 - 


2 /'L +1 f u fh +1 r x C H 

)=- du dv dx _ H dy 

J t 1 ? J u—T/n J tV; J x — T/n V ^ V) 


rj-iZH 

=PhU - k )ff2H- 


where 


Ph(t ) = - ( — \r - 2 \ zh + 4|r — 1| 2H — 6|r| 2H + 4|r + 1| 2H — \r + 2| 2H ), r £ N U {0}. 


Note that (s + t) ^ |s — t\ and 

dv[^d X r 

Jt" Ju-T/n Jt" Jx-T/n \9s 2 8t 2 


I j S H ,n TB H ,n I 

I d jk -d jk ’\ = 


Ch 


(u - yY~ 2H 


dy 


[tj+l 

nu 

Gfc+i 

nx 

dy 

/ du 

/ dv 

/ dx 

/ 

Jtv- 

3 

>u — T/n 

4 

' x — T/n 

(v + y) 4 ~ 2H 

/■*”+! 

nu 

f*k +1 

nx 

dy 

/ du 

/ dv 

/ da; 

/ 

Jtv- 

3 

1 u—T/n 

% 

' x — T/n 

(v — y) 4 ~ 2H 


=| Ch 
<| C H 

p\ PH {j - k)\{Tn~y H 
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for \j — fc| > 1. Thus 


\d^!k ’”| < 2 \pn{j - k)\{Tn x f H for \j - k\ ^ 1. 

It still remains to prove the cases when j = 1 and k > 1 or k = 1 and j > 1. Set ,E = t^—Te/n 
and u e = u + Te/n. Since 


ils" 

rufc 


' — df 


= lim 

e—>0 


<|Ch| lim 

e— >0 


f 2 du [ 

J t™ J U 

f 2 du [ 
JtV 1 J u 


dv 


u £ —T/n 


r t k +1 r x 
/ dx 
J Jx—T/i 

t*I. 


( d 4 G H 
V ds 2 dt 2 


(v,y) + 


dy 


Ch 


(u - yY 


dy 


(v - y) 4 


= hTn- 1 ) 211 lim I [{k + 1) - (2 - e)] 2H - 2 • k 2H + {k + 1 - ef H - 2 [k - (2 - e)} 2H 

2 e—¥ 0 

+ 4 {k - l) 2H - 2 (k - e) 2H + [(k - 1) - (2 - e)] 2H - 2 (k - 2) 2H + (k - 1 - e) 2 " 
=(Tn- 1 ) 2H |p H (fc-l)| 


then the inequality 

|di,k| ^2|p H (fc-l)|(Tn- 1 ) 2H 

holds. A similar argument yields 


| < 2\p H (J ~ l)|(Tn _1 ) 2H . 

Now we shall prove the statement of the lemma. We will use the estimate |8| and the 
equalities Ph{— r) = Ph(v), pn( 0) = 4 — 2 2H , and pn( 1) = —1(7 — 4 • 2 2H + 3 2H ). Note that 
for H ^ 1/2 


max 

1 <Kn 


: V|4"U^ + 2 max V 


j =i 




|p»(j')I _ 7 , 0 7- 4-2 2ff + 3 2 " 


6 ^ I Ph (0)| 6 

In [5] it was proven that 


= v +2 


4 — 4 H 


+*E 


\phU)\ 

4 IM0)l ■ 


E 

3=2 




3 - 3 • 2 2H + 3 2H 
4-4 H 


Thus 


max 

l<fcsCn 


71 1 ,s H 1^7 , n 7-4 - 2 2ff + 3 2H 


uEK 

I=i 


'jk 


^77 + 2 
0 


4-4» 


- 2 ■ sgn(2I7 - 1) 


3 — 3 • 2 2H + 3 2H 
4-4» 


- I 6 


I+ 2 10 - 7 4 4 _4 2 - 32g for H < 1/2, 
1+2 


< 


for 1/2 < H <1 

i, 

1 + 2 for 1/2 < H < 1. 


l + f for H < 1/2, 

7 


3.2 Bifractional Brownian motion 

Definition 10 (pT]) A bifractional Brownian motion (bifBm) B KH = (B^ H ,t ^ 0) with 
parameters H £ (0,1) and K £ (0,1] is a centered Gaussian process with the covariance 
function 

F KH (t,s) = 2~ K {(t 2H +s 2H ) K -\t-s\ 2HK ), s,t> 0. 

The incremental variance function of bifBm is 

2 -A nH,K |2 0 l-ifri, 2 HK /,2 H . 2H\K1 , ,2 HK . 2 HK 

a B H,K{s,t) = tu\B t ’ — B g ’ \ = 2 [|£ — s| — (t +s J J+t + s 
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( 10 ) 


Let H G (0,1) and K G (0,1]. Then 

2 |f — s| ^ ct b h,k (s, t) ^ 2 |t — s| 

for all s,t G [0,oo) (see [TT]). 

It is known that for bifBm the Orey index is equal to HK (see [IB]). If K = 1 then bifBm 
becomes fBm, hence we ignore this case. 

Lemma 11 Assume that B H ' K = {B H,K (t) : t € [0, T]} is a bifBm with K G (0,1) and 
H G (0,1/2). Then 


n — 1 


£4 


- i 




1/2 


22 


where 


ch 


/I n 
\6T n — 1 ' 9In4 

E(A^S ir,ic ) 2 . 


2 l-X T 2if(4 _ 2 2H) 

Remark 12 Without this restriction for H the expressions becomes more complicated. 

Proof. The proof of the lemma follows the outlines of the proof of Lemma [ 8 ] Observe that 
the following equality 


it 

d h 


= 1 - 


b(k, H, K) 
4 - 4 kh 


holds, where 


b(k, H, K) =2 [(k + l) 2H + k 2H ] K + 2 [k 2H + {k- l) 2H ] K 

- 2 K ~ 1 [(k + 1) 2KH + 4 k 2KH + (k - 1) 2KH ) - \{k + 1) 2H + (k- 1 ) 2H Y 


By computer modeling we obtain inequalities 

b kh 

max d kk ^ 1 and 
H,K£( 0,1) 

for all k ^ 1. Let ifG’f and r n = [<p(Tn~ 1 )nT~ 1 ]. Then 


b(k, H, K) 1 

h,S,i) 4 - 2 2KH 6 


( 11 ) 




n — 1 


- 1 


i AbftiU) 
n- 1 ^ 4 -^^ + 


k =1 

=/l+/ 2 . 


b{k,H,K) 

2 2KH ' Tn+1 ^„_! 4 _ 2 2KH 


and 


r ^ T n b(k,H,K) 1 — 1 \rri — 1 « 

/i ^ - max —-^ - tP\ Tn )1 - 

n - 1 i^r n 4 - 6 ' n - 1 

\K \ | iK 


Using the inequality (a + b) ^ 2 K ( a K + 6 K ), a, b ^ 0, 0 < K < 1, and the formula 


(i+*)“ = i+ 

k =1 

we obtain 
| 6 (fc, H, A')| =k 2HK 


a{a — 1 ) • • • (a — k + 1 ) k 
k\ * ’ 


1 + l 


+ 2 




<2K 


- 2 


=2k 


- 2 


2HK I 0 1-Jf 


1 + 


+ 1 


— 2 


1+ i 


for — 1 < x < 1 , 


+4+(1 -fc 


1+ ^ 


2 H 


1 + l 


k 

2 KH 


+ 2 +| 1 - 

1 


1 _ fe 


+ 2+(1 -fc 


k 

2 KH 


1 2 1_x 

', , ^ 2H(2H - 1) (2H - 2m + 1) ( 

[ \ 2 m 

r 

^ ( 2 m)! V 

L m= 1 V y 

fc/ 


4+E 


2KH(2KH -!)••■ {2KH - 2m + 1) / 1 \ 2m 


(2m)! 


(O' 


:= b(k,H,K). 
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Note that the sign of 2KH(2KH — 1)(2KH — 2) • ■ • ( 2KH — 2m + 1) is the same as that of 
2KH — 1. Thus the estimate of b(k, H, K) depends on the signs of 2H — 1 and 2 KH — 1. Then 
for H G (0,1/2) and k ^ 2 


Thus 


b{k , H, K) ^ 2 K k 2HK ^ 2KH{1 - 2I\H ) • • ■ (2m - 1 - 2ATf) / 1 \ 2 ™ 
4-4 KH ^4-4 kh ^ (2m)\ \k) 

m= 1 v 

< 2 K k 2HK {l- KH) y, j_ / 1 N2m 2 K k 2HK y /l\ 2m 
4 _ m \k) ' In4 m \k) 

m =1 m=l 

2 K k 2HK / 1 \ 11 • 2 k 1 22 x 

In 4 k 2 )^ 9 In 4 k 2 ~ 2KH ^ 9 In 4 



E* 





n 22 T 
n — 1 91n4 <p(Tn~ 1 )n 


Let <p(Tn x ) = (Tn 1 ) 1 ' 2 . It belongs to the class of functions 4b After plugging it into the 
obtained inequality we get 

-j—vV™- 1 A( T \ 1/2 fi 22 T 

n — 1 jy kk ^ 6 \n) T(n — 1) 9 In 4 (Tn~ 1 ) 1 / 2 n 

<1/J_y /2 _n_ 22 /T\ 1/2 _ /T\ 1/2 / 1 n 22 \ 

^6 \ nT) n — l^~91n4\n/ \ n / \6 T n — 1 + 9 In 4 / 


Lemma 13 Assume that B H ' K = {B H,K (t) : t G [0, T]} is a bifBm with K G (0,1) and 
H G (0,1/2). Then 

n — 1 

1 <k<n— 1 51 \ d f* H I < 5 - 005 - ( 12 ) 

" " n j =1 

Proof. The fourth order mixed partial derivative of the covariance function R KH (s,t) has 
the form 


8 4 Fkh 

ds 2 dt 2 


(M) 


C' 


(i) 


1 ( 2 ) 


|S — t|4-2 KH + C KH 


, ,\4if —2 ( 2 H . ,2H\K-4 

(St) (s +t ) 


+ C^ H (st) 2H ~ 2 (s 2H +t 2H Y 


for each s, t > 0 such that s t, where 

C^h —2HK(2KH - 1)(2 HK - 2)(2 HK - 3)2 _/f , 
Cgh =K(I< - 1 ){K - 2){K - 3)(2H) 4 2~ K , 

C { kh =K(I< - l){2Hf{2H - 1)2 ~ k (2KH -2 H - 1). 


Since 2 s H t H ^ s 2H + t 2H and K — 2 < 0, K — 4 < 0 it follows that 


(. st) 2H ~ 2 (s 2H +t 2H ) K ~ 2 ^2 K - 2 (st) KH ~ 2 , 
(st) 4H ~ 2 (s 2H + t 2H ) K ~ 4 <2 K ~ 4 (st) KH ~ 2 . 


(13) 

(14) 


Let B kh be a fBm with the Orey index KH. Assume that KH ^ 1/2 and j / 1 or 1; / 1. 
Using the inequalities (131, (14) and 

T 

vy=[y + (v- y)]y > [v - y)y > (v - y) -, y ^ t", 


10 
































we obtain 


■jK H 




I df k — 2 dj k 
< ir (2) 

*5 \ L KH 


+ \C\ 


f l j +1 f u rk+i r 

/ du dv dx 

J t™ Ju — T/n Jt 7 ^ Jx 

/•*?+! r u /■*£+1 r 

/ du dv dx 

Jt 1 } Ju — T/n J t 7 ^ Jx 

tu 


dy 


—T/n {vy) 2 ~^{v 2H +y 2H Y 
dy 


nt+i f u f f h+i r 

/ du dv dx 

J t 7 } Ju — T/n Jt™ Jx 

t 


_r/n {vy) 2 ~ 2H {v 2H + y 2H ) 2 ~ K 

2 K - 2 \C ( / ) H \+2 K - A \d ( / ) h 




x — T/n 


( vy) 2 ~ KH ' dy 

2 /f - 2 |C^|+2 A '- 4 | C$ H 

/ x — T/n 


(v - y) 2 


f f j +1 /‘ t fc + l f 

/ du dv dx 

Jt 1 } Ju — T/n 

2—KH r\K—2 1 ^(2) | . r\K— 41 ^(3) i 

2 |g ""' + 2 , |p™(* - fc)| {Tn~ 1 ) KH+2 

\^KH\ 

2KH r\K — 2 | s-i(2) | , c^K— 4|p}(3) i 

2 ICWI + 2 l“^£i|p Jfff( j-fc )| , 


dy 


\Ckh\ 

for | j — fc| ^ 1, where 

Ckh =(KH - 1 )KH(KH + 1)(KH + 2) 

P*h(? - fc) =6(j - k ) KH+2 - 40 - k - 1) KH+2 - 4(j - k + l) KH+2 
+ (j -k- 2) kh+2 + (j — k + 2) kh+2 . 

So 


I ^B kh I 

I dj k 


^2^\p KH (j-k)\[/) + 


T 


2KH 2 K - 2 \C( 2) H \+2 K - A \df ) 


I c* 


Ipk-hC? - fe )l ( - 

n 


The cases when j = 1 and A; > 1 or k = 1 and j > 1 can be proven in a way analogous to that 
of subfBm. Next, we obtain 

n—1 n—1 

, max V |df fc | sj 1 + ( max |df+" fc | + max |d*_ ljJfe |) + max V | df k \ . 


2 = 1 


2 = 1 

|j —fc|S*2 


Since 


\dC H k \ Jd£ H dg* k+1 , 


then by the inequality (111 we get 


b kh 

max |d fc+ i fc | ^ 1. 

l^.k^.n — 2 


Reasoning as in [5j, Appendix A we obtain 


\PKH(j)\ = 3 - 3 • 2 kh+2 + 3 kh+2 . 

3=2 


Furthermore, 


max 


L, E K 


'jk 


3 = 1 

|j-fc|^2 


E 


IpxhQ' - fc)| 
4 _ 41 m 


< max £ l T ( i^ )l + 23iC " 3|eg ^ +2 r 5|ag ”' — 

l^fc^n-1 ' 4_4«n \C KH \ l<fc<n-l 

J=1 11 J=1 

li-fc|>2 |j-fc|>2 

1 1 , 3 - 3 • 2 2KH + 3 2KH 2 2K - 3 \d i / ) H \ + 2 2K - 5 \C ( / ) H \ 3-3-2 kh+2 + 3 kh+2 

<-- sgn(2A_ff - 1)-^-+-j^j-- 
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It is clear that 

2 2 k- 3 \Q(2) h \ ^ _ K ^ 2 _ K ^ 3 _ k)2 4 H 3 2 k ~ 3 _ (H - HK)(2 - K )(3 - I<)2 4 H 2 2 k ~ 3 


| C* 


and 


(1 - KH)(KH + 1)(KH + 2) 

^ (2 - A") (3 - K)H 2 2 k+1 ^ 6 ■ 2 if_1 ^ 

^ (A'A + 1)(KH + 2) " 2 " 

2 2K- 5 |g(3)^| 2 k_3 (1 - A)A|2tf — 1|(1 + 2H — 2KH) 


(1 - KH)(KH + 1)(KH + 2) 


| G h 


(1 — KH)(KH + 1)(KH + 2) 

^2 /f - 3 |2A- 1|(1 +2H- 2KH) ^ 2-2 K ~ 3 ^ 1 
" (KH + 1)(I<H + 2) " 2 ^ 4 ' 


Thus 


max 

l</c<n- 


i=i 


'jk 


<3 + 


3 _ 3. 2 2ifH + 3 2KH 13 3 - 3 • 2 X "+ 2 + 3 KH + 2 


2(4 — A KH ) 


+ T 


4 _ 4 kh 


< 3+ j + t® = 5 ' 005 


since the numerator of the second term is a decreasing function, and the numerator of the 
third term is an increasing function. 

3.3 Ornstein-Uhlenbeck process 

The fractional Ornstein-Uhlenbeck (fO-U) process of the first kind is the unique solution 
of the stochastic differential equation 


J" 

Jo 


X t = xo — n / X a ds + 9B, 


t € T, 


(15) 


with fi, 9 > 0, where B H , 0 < H < 1, is a fBm. Its explicit solution is given by 


X t = x 0 e _Mt +9 [ 

Jo 


e-^-^dB”, 


where the integral exists as a Riemann-Stieltjes integral for all t > 0 (see, e.g., [6]). First we 
show the following lemma. 

Lemma 14 Let X be the solution of equation 15). Assume that B H = {B H (t) : t £ [0,T]} 
is a fBm with H £ (0, H*], where a real number H* < 1 is known. Then for H ^ 1/2 


max 

1 


LjZK*’" 




4 iiT 


3 = 1 


4 - 2 2ff * 


^(2p- 


Cl* 


+1 


+ V2(4-2 2 »*) 




+1 


1/2 


where C = p 2 [3 Xq9 2 + 6T 2 ] and 

„2 H 


lX.n _ n J. X,n 

a jk ~ Q2J'2H ^4 _ 2 2h ) a J k ’ 


df k n = E(A^XA^X). 


n,j n,k y 


A n , k Z — fi j X s 


ds. 


Proof. Denote 

It is clear that 

\d*k n - 9 2 df k H ’ n | =\E[A n , j+1 Z - A B , J -Z][A„, fc + 1 Z - A n , k Z\ + 9E[A nJ+1 Z - A nJ Z]A™ k B H 
+ 9E[A n , k+1 Z - A n , k Z]A^B H \ 

and 


[A n j + iZ A U) jZ] — p 


fj+i r l j 

J [X s - Xj] ds- J [X s - Xj] ds 
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Reasoning as in m we obtain 


supE.Y 2 «: 3x 2 0 + 69 2 T 2 


and 


ElXt-Xj] 2 s?2/i 2 ^^ EX 2 ds + 20 2 (J^j ^2n 2 [3xl + 69 2 T 2 ](J^j + 29' 


(T^ 

2 H 

t \ 2 ~- h 

- 


c - +1 

\nj 


L \ n ) J 


where C = /x 2 \?>Xq6 2 + 6T 2 ]. Thus 


I n /q2 ,7i I o 2 /)2 ( 

I d 7fc d 7fc 61 ( “ 


2+2H 


c| s' +I 


+ 4^2(4-2 2 «)/i6) 2 ^ 


l+2iT 


c 


2-2 H -| 1/2 

+ 1 


Consequently, 






max 

J = 1 

8/1 2 T 2 
: 4- 2 2H " "rT 

4 /iT 


+ max 

l<fc< 7 i 


2 ,71 


n 


4^2 /iT 


j'=i 


71 — 1 

\ ' I jB H ,77 

+ max > d. fc 
1 = 1 


+ 1 


n 


( T 

„ T\ 

2/1 — 

C - +1 

V n 

L \ n ) \ 


V4 - 2 2 » 

+ \/2(4 — 2 2H *) 


2-2 H* 


+ 1 


1/2 


+ 3 


4- 2 2 » 

Lemma 15 Let X be the solution of the equation Il5l). TTien 


C(I 

n 


+ 1 


1/2 


+ 3‘ 


1 71—1 

*- jX,ti -> 

_ 1 / j ^kk 1 


n — 1 


€ 


4 - 2 2 »* 

+ V2(4-2 2 »’) 


T 

( T 


/ rjn \ 2 — 2 H 

- 

2/i — 

c 

— ) +1 

n 

V n 


K n J J 


Proof. To prove this lemma, observe that 

71—1 

n — 1 


1 n_1 
1 rX.Tl 

_ 1 / ^ ^ 


^ I iX.n iB r 

^I™<n-1 ^ kk ~ dkk 


C T 

n 


+ 


+ 1 


1/2 


i n_1 

1 ,b h 

—7 ® fcfe 


n — 1 


- 1 


I iX ,n -t I 

= max \duu — 1 . 

l<fc<n-l 1 1 


From the inequality 


I ^ /}2 7B n i'll I q 2/)2 / 

“fcfc ~ 0 d <8/x 0 — 


2+214 


n 


+ 1 


/ 

1+214 

/ t\ 2 ~ 2H 

- 


c - +1 

Vn, 


L \ n J \ 


1/2 


it follows that the statement of the lemma holds. 

4 Simulations 

The simulations of the obtained confidence intervals presented below were performed using the 
R software environment [16]. Sample paths of fBm were generated using the circulant matrix 
embedding method and were further used to simulate the sample paths of the fractional 


Ornstein-Uhlenbeck process (151. The constants for the latter were (arbitrarily) chosen as 
xo = 0 and /i = 0.5. Sample paths of the sub-fractional and bifractional Brownian motion 
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were simulated using the Cholesky method. Due to the notable computational requirements 
of this method the maximum sample path length considered was n = 1600. Figures presented 
below correspond to the case of the confidence level 1 — a, a = 0.1. The observed coverage 
percentages in all cases were at least as good as claimed in Theorem [6] 

Figures[l]-[5]present the confidence interval (Cl) lengths for all the process types considered 
in this paper. Figure [6] shows the median ratios of the confidence intervals lengths, where the 
Cl lengths of the subfBm, bifBm and fO-U processes were divided by the corresponding Cl 
lengths of fBm. It can be seen that in almost all cases the confidence intervals behave in a 
similar way, one notable exception being the case of fO-U as the value of H approaches 1. 
This is hardly unexpected given the normalization used in Lemma [14] and in this scenario the 
Cl covers the whole interval of possible parameter values 0 < H < 1. 


fBm Cl length (alpha=0.1) 



Figure 1: I H%»{a) - H l n n f (a) I 


subfBm Cl length (alpha=0.1) 



Figure 2: I H™P(a) - H l n n f(a) I 


bifBm Cl length (alpha=0.1 ,K=0.1) 


bifBm Cl length (alpha=0.1 ,K=0.9) 




n 


n 


Figure 3: \(HK) s n u P(a) - 


Figure 4: I (HK) s n u P(a) - {HK)™f(a) 
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0-U Cl length (alpha=0.1) 


Ratios of Cl lengths 



n 


— bifBm / fBm 

- - subfBm / fBm 

O-U/fBm 



200 400 600 800 1000 1400 

Sample path length 


Figure 5: | H s n u P{a) - 


Figure 6: Comparison of Cl lengths 
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Appendix. Code listings 

genFBM.r 

# genFBM returns a single sample path of the fBm # 

# H — the Hurst index # 

# N — length of the sample path # 

genFBM <— function (H,N) { 

#M— length of the required vector of autocovariances 
M <— 2'(ceiling( 1 + log (N—l)/log (2) )) 

# V is the M-vector of autocovariances of fGn 
H2 <- 2*H 

V <- c(0:(M/2) , (M/2 — 1): 1) 

V <- l/(2*N~(H2))*(abs(V-l)~(H2) - 2*V~(H2) + (V+1)~(H2)) 

# W is the fast Fourier transform of V 
W<- Re( fft (V)) 

# We increase M until all coordinates of W are positive 
while (any(W<=0) & (M<2'30)) { 

M <- 2*M 

V <- c (0: (M/2) , (M/ 2 — 1): 1) 

V <- 1 / (2 *N" (H2)) * (abs (V— 1) * (H2) - 2*V" (H2) + (V+1)~(H2)) 
W<- Re( fft (V)) 

} 

#X and Y and iid Gaussian with mean 0 and variance l/sqrt(2) 
X <— rnorm(M, mean=0, sd = (l/sqrt (2))) 

V <— rnorm(M, mean=0, sd = (l/sqrt (2))) 

Z <— vector (length=M) 

Z[l] <- X[l] 

Z [M/2 + 1] <— Y[l] 

Z [c ( 2 : (M/ 2 )) ] <- X[c(2:(M/2))] + 1 i *Y[c (2 : (M/2))] 

Z [M+2—c ( 2 : (M/ 2 )) ] <— X [ c (2: (M/ 2)) ] - 1 i *Y[ c (2 : (M/ 2)) ] 

U <— sqrt(W)*Z 

# X is the fast Fourier transform of U = sqrt(W) * ( X + iY ) 

X <— fft(U, inverse=T) 

# The real part of the first N coordinates of X is the desired 

# sample path, in this case, the fGn. 

# We calculate the fBm sample path as cumulated sums of fGn. 
BHinc <- c(0,Re(X[l:N]))/sqrt(M) 

BH <— cumsum( BHinc) 
return (BH[ — length (BH) ] ) 

} 

Cholesky.r 

# genSFBMch generates a batch of sample paths # 

# of sub-fractional Brownian motion # 

# H — the Hurst index # 

# N — length of sample paths # 
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# Nsp — number of sample paths # 


genSFBMcli <— function (H,N, Nsp) { 

# sfBm’s covariation matrix 

sfBm_cov <— function (t , s ,H) { 

H2 <- H*2 

cts <- s"(H2)+t“(H2)-((t+s)*(H2)+abs(t-s)'(H2))/2 
return( cts ) 

} 

G <— matrix ( ncol=N, nrow=N) 
for (j in c(1:N)) { 
i <— seq (1 :N) /N 

G[,j] <— (apply (as . matrix (i ), 1 , sfBm _cov , j/N,H)) 

} 

# Cholesky decomposition 
L <- t ( chol (G)) 

# a function to generate a single sample path 
gen_sp <— function(L) { 

Z <— rnorm(N) 

B <— c (0 , (L %*% Z)) 

return(B) 

} 

# generating Nsp of them 

BHM <— replicate(Nsp,gen_sp (L)) 
return (BHM) 

} 

mmMmntmffffm/f/mMmmmmMmm f# 

# genSFBMch generates a batch of sample paths # 


# of bifractional Brownian motion # 

# H, K — parameters of the bifBm # 

# N — length of sample paths # 

# Nsp — number of sample paths # 


genBFBMch <— function (H,K,N, Nsp) { 

# bifBm ’ s covariation matrix 

bifBrn_cov <— function (t , s ,H) { 

H2 <- H*2 

cts <- ((abs (t)' (H2) + abs ( s ) ~ (H2))' K—abs (t—s ) * (K*H2)) / (2 'K) 
return ( cts ) 

} 

G <— matrix ( ncol=N, nrow=N) 
for (j in c(1:N)) { 
i <— seq (1 :N) /N 

G[ , j ] <— (apply (as . matrix ( i ), 1 , bifBm _cov , j/N,H)) 

} 

# Cholesky decomposition 
L <- t(chol (G) ) 

# a function to generate a single sample path 
gen_sp <— function(L) { 

Z <— rnorm(N) 

B <— c (0 , (L %*% Z)) 
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return(B) 


} 

# generating Nsp of them 

BHM <— replicate(Nsp,gen_sp(L)) 

return (BHM) 

} 

fbmCI.r 

The following code evaluates the confidence intervals for the fractional Brownian motion. 
Simulations for other processes considered in this paper were performed in a similar way. 

source (” Cholesky . r” ) 
source (’’genFBM. r” ) 

mmsmMmmmmmm 

# D efinitions of functions # 

mmsmMmmmmmm 

phiL = function(z, nu, eps , N) { 

if (z >= 0 &: z < (eps + 1) * sqrt (N— 1)) { 
return ( exp(z*sqrt (N—1) / (2*nu)) 

* (1 — z/((eps + l) * sqrt (N— 1))) 

**( (eps + 1) * (N—1) / (2*nu) ) ) } 
else return(O) 

} 

phiR = function(z, nu, eps, N) { 
return ( exp(—z*sqrt (N—1)/ (2*nu)) 

* (1 + z/((eps + l) * sqrt (N— 1))) 

**( (eps + 1) * (N—1) / (2*nu) ) ) 

} 

pliiLOpt = function (z, nu, eps, N, value){ 
return ( abs ( phiL ( z , nu , eps, N) — value)) 

} 

phiROpt = function (z, nu, eps, N, value){ 
return ( abs (phiR(z, nu, eps, N) — value)) 

} 

gn = function (x, N) { 

if (x = 0) { return(— Inf) } else { return( 2 *x*log (N) — log(4 — 4**x)) } 

} 

gnOpt = function (x, N, value) { 
return(abs( gn (x , N) — value)) 

} 

logStar = function(x) { 

if (x>0) { return ( log(x )) } else { return(— Inf) } 

} 


confBC = function(B, alpha, nu, eps, N, kappa) { 

qL = optimize ( f=phiLOpt , interval=c (0 , 20), value=alpha/2 , 
nu=nu , eps=eps , N=(N— 1)) [ [ 1 ] ] 

qR = optimize ( f=phiROpt , interval=c (0 , 20), value=alpha/2 , 
nu=nu , eps=eps , N=(N—1 )) [ [ 1 ] ] 

xL = 1—qL/sqrt (N—1) 
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xR = 1+qR/sqrt (N— 1) 

Sn = mean( diff (B, differences =2)**2) 

argL = max( logStar((xL — eps)*kappa**2 / Sn), —log(3) ) 
argR = log((xR + eps)*kappa**2 / Sn) 

Hmin = max( c (0 , optimize (gnOpt , c(0.00001, 0.99999), 

N=N, value=argL ) [ [1]])) 

Hmax = min( c (1 , optimize (gnOpt, c(0.00001, 0.99999), 

N=N, value=argR)[[l]])) 

return (c (Hmin, Hmax)) 

} 

# Parameters used for modelling # 

nseq = c(200, 400, 800, 1600) 

Hseq = c (0.01, 0.25, 0.75, 0.99) 

Nsp = 1000 
alpha = 0.1 

kappa = 1 

# ciLenM — matrix of Cl lengths # 

# cilnsideM — matrix of Cl coverage ratios # 

mMsmmsmmmsmsmmmMmsmm 

ciLenM = matrix (nrow=length( Hseq) , ncol=length ( nseq )) 
cilnsideM = matrix (nrow=length (Hseq) , ncol=length ( nseq )) 
for (m in 1: length (Hseq)) { 

H = Hseq [m] 

ciLen <— vector (length=length ( nseq )) 
cilnside <— vector (length=length ( nseq )) 
for (k in 1:length( nseq )) { 

N = nseq[k] 

BM = replicate (Nsp, genFBM(H, N)) 
eps = 1/N 
nu = 8/3 

Cl = apply (BM, 2, confBC , alpha=alpha , nu=nu, eps=eps , 
N=N, kappa—kappa ) 
ciLen [k] <— mean(CI [2 ,] — Cl [ 1 ,] ) 

cilnside [k] <— length (which (Cl [ 1 ,] <H & CI[2,]>H))/Nsp 

} 

ciLenM [m,] = round( ciLen, 3) 
cilnsideM [m, ] = cilnside 


# Plotting the results # 

//////////////////////////////////////////m 

pdf (”bw_fbm_ci . pdf” , width = 8, height = 4) 
par (mfrow=c (1 ,2)) 
yl = c ( min (ciLenM ) , max(c.iLenM)) 

plot(nseq, ciLenM [1,], type=”b”, lwd = 2, lty=l, xlab=”n” , ylab=””, 

main=paste0 (”fBm„CI„ lengt h „( alpha—’ , alpha,”)”), ylim=yl , xaxt=”n”) 
for (k in 2 : length (Hseq)) { 
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lines (nseq, ciLenM[k,] , type=”b” , lwd = 2, lty=k) 

} 

abline(h=0, lty=2) 

legend (” topright ” , c(”H=0.01”, ”H=0.25” , ”H=0.75” , 
lwd = 2, 1 ty =1: length ( Hseq )) 
axis(l, nseq, labels=nseq , las=2) 

yl = c (min( cilnsideM ) , max( cilnsideM )) 
plot(nseq, cilnsideM [ 1 ,] , type—’b” , 11y = 1, xlab=”n” 
main—’ Coverage„%” , ylim=yl , xaxt-’n” ) 
for (k in 2 : length (Hseq)) { 

lines (nseq, cilnsideM [ k , ] , type=”b” , lty=k) 

} 

legend (” bottomright ” ,c(”H=0.01” , ”H=0.25” , ”H=0.75” 
lwd = 2, 1 ty =1: length (Hseq)) 
axis(l, nseq, labels=nseq , las=2) 
dev . off () 


H=0.99” ) , 


, ylab—’ ” , 


, ”H=0.99” ) , 
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